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Dedicated to Percy Deift with gratitude and admiration. 

Abstract. In this paper we derive uniform asymptotic expansions for the 
partial sums of the exponential series. We indicate how this information will be 
used in a later publication to obtain full and explicitly computable asymptotic 
expansions with error bounds for all zeros of the Taylor polynomials pn-i{z) = 
S^Zq z^/ fc! ■ Our proof is based on a representation of p„_i(nz) in terms of 

an integral of the form f "^^ ^ ds. We demonstrate how to derive uniform 
expansions for such integrals using a Riemann-Hilbert approach. A comparison 
with classical steepest descent analysis shows the advantages of the Riemann- 
Hilbert analysis in particular for points 2 that are close to the critical points 
of <b. 



1. Introduction 

During the past fifteen years and largely due to the ground breaking work [5], 
[5] of Deift and Zhou, Riemann-Hilbert problems have become a powerful tool in 
asymptotic analysis with applications in many fields such as inverse scattering the- 
ory, integrable PDE's, orthogonal polynomials, statistical mechanics and random 
matrix theory. With this paper we would like to add another item to this list of 
applications, namely the study of the asymptotic behavior of zeros of Taylor poly- 
nomials of entire functions. Although our method works in principle for a large 
class of entire functions we will restrict our attention to the classic case of Taylor 
polynomials of in order to keep the presentation as non technical as possible. 
It is also one of the intentions of this paper to highlight the advantageous aspects 
of the Riemann-Hilbert approach in a simple situation. The first and maybe most 
surprising feature is how many objects of interest can be characterized as the unique 
solution of some Riemann-Hilbert problem and with this paper we present one more 
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Figure 1 . Szego curve (solid line) and rescaled zeros of p2o (tri- 
angles) and pso (boxes) 

illustration of this fact. Secondly, there is a set of techniques to transform Riemann- 
Hilbert problems to some form where an asymptotic expansion of the quantity of 
interest can be extracted in a systematic fashion. In our case, we will only need 
one of these techniques and that is the construction of a local parametrix. This 
construction exemplifies in a very simple situation just one of the wonderful ideas 
of Deift and Zhou [5], [6], namely how to effectively localize a Riemann-Hilbert 
problem which by definition is a global problem in the complex plane. Localization 
is an indispensable step in the asymptotic analysis of any Riemann-Hilbert problem 
where the asymptotic behavior of the solution is different for different regions of 
the complex plane. Moreover, and that is the third aspect we want to highlight, 
the localization procedure of Deift and Zhou always provides for matching asymp- 
totics in those regions where the asymptotic behavior changes and one obtains error 
bounds uniform in all of C. 

We denote by p„(z) := 1 + z + ■ ■ ■ + ^ the partial sums of the exponential 
series. The problem to describe the asymptotic distribution of the zeros of Pn 
was posed and solved in the classical paper of Szego [TT]. He proved that the 
zeros of Pm divided by n, converge in the limit n — > oo to some curve -Doo, now 
called Szego curve, which consists of all complex numbers l^l < 1 that satisfy the 
equation \ze^~^\ = 1. The curve Doo together with the rescaled zeros of pn for 
n — 20 and n — 80 are displayed in Figure [TJ Moreover, Szego also determined the 
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limiting distribution of the rescaled zeros on Dqo- A first result presenting error 
bounds on the distance between the Szego curve and the zeros of p„(nz) has been 
established by Buckholtz ^ who showed that they are located in the exterior of 
Z?oo at a distance of at most 2e/^/n. For each zero of its distance from Doo is 

measured by the minimal distance between the zero and all the points of the Szego 
curve. Subsequently more detailed asymptotics of pn{nz) have been derived. It 
turns out that zq = 1 is a critical point where the asymptotic behavior changes. It 
was shown by Newman and Rivlin [9J that in neighborhoods of zq of size 0(1 /y/n) 
the rescaled polynomials p„ (n + Wy/n) can asymptotically be expressed in terms of 
the complementary error function erfc. Carpenter, Varga, and Waldvogel ^ (see 
also [13] and see [12] for an interesting discussion of zeros of the partial sums of e^) 
then provided an asymptotic expansion oi pn{nz) in compact subsets of C \ {zq\. 
These results were used in f4j to obtain lower and upper bounds on the distance of 
the zeros from the Szego curve. Note that up to the recent paper by Bleher and 
Mallison 2J no uniform asymptotics for p„(nz) and its zeros were available in a 
fixed size neighborhood of the critical point zq. In a different direction and using 
methods from logarithmic potential theory Andrievskii, Carpenter and Varga [1] 
recently extended the results of Szego ^llj on the angular distribution of the zeros 
by proving uniform error bounds for regions including zq. We refer the reader to 
the recent review [14] for a description of more results on the zeros of p„ . Related 
results on the zeros of the partial sums of cos, sin and of more general sums of 
exponential functions can be found in 2J and in references therein. Finally, we 
mention the review [1^ on the behavior of the zeros of more general sections and 
tails of power series. 

The main novel result proved in the present paper is Theorem l3.3l which provides 
an explicitly computable asymptotic expansion near zq for a quantity Fn{z) closely 
related to Pn-iinz) (see (|1.1|) . (|1.2|) below). Together with Remark l3.4l we produce 
an asymptotic expansion of Fn with uniform error bounds for z contained in some 
region V where all the roots of of j3„_i(nz) are located. These results will be used in 
a subsequent publication to derive explicitly computable asymptotic expansions for 
all the zeros of Pn-i in terms of the zeros of the complementary error function for 
zeros close to the critical point, otherwise in terms of the solutions of (ze^^^)" = 1 
which lie on the Szego curve. We will state the corresponding results without proof 
in Theorems 14.11 and 14.21 below. 

Let 7 be any smooth Jordan curve encircling the origin counterclockwise. For 
z e C \ 7 define 

1 f s~"e"^*'^"^' 1 f e"'^''*-' 

(1.1) Fjz) := / ds = / ds, with 0(s) = s - 1 - Ins, 

zni s — z 27ri s — z 

where In denotes the standard branch of the logarithm slit along the negative real 
axis. A straightforward application of the calculus of residues then yields 

, . , , „ , , _ J — F„(z) , for z in the exterior of 7 

(1.2) (ez) p„_i(nz) = j ^„^(^) _ ^^^^^ ^ for z 7^ in the interior of 7 

To the best of our knowledge the Cauchy-type representation of p„_i provided by 
(|l.ip and (|1.2p which is crucial for our Riemann-Hilbert approach has so far not 
been used in the asymptotic analysis of the zeros of 

Given the form of Fn it is natural to use the method of steepest descent for 
the asymptotic analysis of the integral. This method requires that the path 7 of 
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integration passes through the critical point zq = 1 of the function (p. This indicates 
already that the asymptotic analysis will be most difficult for z close to 1 due to 
the term (s — z)~^ in the integrand. We discuss this approach in detail in Section 
[2] below and we will see that by standard techniques of steepest descent analysis 
one may derive an asymptotic expansion for Fn{z) with uniform error bounds for 
|z- 1| > n"" and < a < 1/2. 

Section [3] contains the central result of this paper where we derive an asymptotic 
expansion for Fn in some fixed neighborhood C/e(l) of the critical point. Here we 
use that Fn can easily be characterized as the unique solution of a Riemann-Hilbert 
problem since by definition F„ is just the Cauchy transform of e"*^ (see {RHP)i 
described at the beginning of Section [3]). Using the standard Riemann-Hilbert 
technique of constructing a local parametrix near the critical point the derivation 
of the asymptotic expansion is no more difficult than the computations in the non 
critical situation of Section [21 

In the final Section [4] we briefly explain how to use the asymptotic results to 
obtain asymptotic information on the location of the zeros. This, however, is a 
somewhat technical affair and we will present a complete account of all results and 
their proofs in a later publication. In Theorems 14.111421 we state our results for the 
zeros of Pn-i in the upper half plane. 



In order to remind ourselves of the method of steepest descent (see also [8] for 
an elementary exposition) we first determine the large n asymptotics of a quantity 
that is related to F„ defined in p.ip but is somewhat simpler to analyze. Let 



where 7 denotes some smooth Jordan curve that is oriented counterclockwise con- 
taining the origin in its interior. 

The method of steepest descent offers a recipe how to deform the contour of 
integration in such a way that the asymptotic behavior of G„ can be determined 
most conveniently. More precisely, we are to consider contours that pass through 
critical points of (j) along the path of steepest descent with respect to the real part of 
4>. Such a choice of the contour ensures that the modulus of the integrand obtains its 
local maximums only at the critical points of (j) and - up to exponentially small error 
terms - the integral is determined by the contributions of those parts of the contour 
which lie in small neighborhoods of the critical points. In the situation at hand the 
function (j) has only one critical point at zq = 1 and we can easily understand the 
behavior of the real part of (j). It has a saddle point at zo = 1 and the solid curve 
in Figure m displays those numbers z satisfying 3fJ(0(z)) = 3fJ((/)(zo)) = 0. Note, 
that the Szego curve D^o which was displayed in Figure [H above coincides with the 
closed loop part of the solid curve. Moreover, the dotted line in this figure shows the 
path of steepest descent away from the saddle. This path can easily be determined 
to consist of those points in C \ {0} satisfying argz = 5z by using the property 
that the imaginary part of <j) remains constant on it, i.e. 3(0(z)) = 9(0(zo)) = 0- 
Choosing the smooth closed contour 7 to coincide with the path of steepest descent 
in the right half plane and with the dashed line in the left half plane (in the left 



2. Classical steepest descent analysis 



(2.1) 
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Figure 2. 



half plane any smooth curve that does not intersect the solid line will do) we obtain 

(2.2) G„ = / e'"f''-'Us = 7^ [ e'"*'^'^ds + O (e-"=) , 

where U is any fixed neighborhood of = f . It is only the number c > in the 
error term that will depend on the choice of U . 

In order to further analyze the integral along 7 n ?7 we change variables near the 
critical point of Since 4>"{zq) ^ we can find an open neighborhood Uq of zq 
and a biholomorphic map A : [/^^(O) — > C/q for some (5o > such that 

(2.3) ^{m)=e iorieUsM 

Observe that A maps the imaginary axis onto the path of steepest descent, more 
precisely A(i]R n [/^^(O)) = 7 fl Uq. This can be seen from the characterization of 
the path of steepest descent by the imaginary part of i.e. 

7 n J7o = e C/o : 3?0(z) = 3f/'(l) = and 5R0(z) < ^<j){l) = 0} . 

Choosing U :— Uq in ()2.2p we obtain 

1 pi^o ^ r^°^ 2 

G„ = — / e"« X'im + (e-"=) - — ^ / e"* X'{it/V^)dt+0 (e""^) 

Expanding A at ^ = leads to an asymptotic expansion for Gn- For example, 
we learn from 

that 

Observe that the calculus of residues implies Gn = e~"n"~"'^/r(n) directly from 
p.ip and we have thus found some version of Stirling's formula. 
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It is clear that we may extend our reasoning to integrals of the form 

(2.4) — / /i(s)e"'^(")ds 

27ri ,L 



ft.(A(^^/^A^))e"* X' {it/y/^)dt. 



for analytic functions h. Up to exponentially small error terms (|2.4p is given by 
1 

This expression can easily be expanded in powers of 1/n using the expansions of 
A around and of h around 1. For example, in the case h{l) ^ we immediately 
obtain 

(2.5) / /i(s)e"^(^)ds = (1 + 0(71-1)) 



for any Jordan curve winding counterclockwise around the origin. 

We are now in a position to apply our reasoning to F„ by replacing the analytic 
function h in (|2.5p by the meromorphic functions hz{s) — {s — z)~^. Indeed, we 
have 

Before formulating a result on the large n behavior of Fn{z) we specify the type of 
contours 7 to be considered. 

Definition 2.1. A contour 7 is said to be admissible if 

(i) 7 is a smooth Jordan curve winding counterclockwise around the origin, 
and 

(ii) 7 has a positive distance from the solid curve in Figure [21 except for a part 
that lies in some neighborhood C/ of zq = 1. In this set U the contour 7 
coincides with the path of steepest descent (dotted line in Figure [2|). 

The steepest descent analysis described above then gives 

Theorem 2.2. For any admissible contour 7 and 2: € C \ 7 we have 

F^{z) = -^=J- (1 + O(n-i)) , 

V27rn(l — z) 

where the error term is uniform for z G C \ (C/c(l) U 7) for any e > 0. 

Proof. It remains to discuss the uniformity of the error bound. The only step in the 
derivation of the asymptotic formula for which the uniformity appears problematic 
is the estimate used in ()2.2|) for points z that lie arbitrarily close to 7 because of the 
(s — z)~^ term in the integrand. However, since z stays away some fixed distance 
e from the critical point zq, one may deform the contour of integration for each z 
in such a way that the value of Fn does not change and such that the deformed 
contour is separated from both, the solid line in Figure [2] and from z by a minimal 
distance that only depends on e. □ 

A somewhat more refined analysis provides uniform error bounds also outside 
shrinking disks |z — 1| > n^" for any fixed < a < 1/2. In this case the error term 
O(n-i) in Theorem 12. 21 has to be replaced by 0{n~'^\l — z|^^). For values of z that 
lie in shrinking discs |z — 1| = 0{n~°') it is still possible but cumbersome to use the 
method described above for computing the asymptotics of Fn{z). It is precisely this 
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situation in which wc want to demonstrate how ideas from the asymptotic analysis 
of Riemann-Hilbert problems can be used. 

3. RiEMANN-HlLBERT ANALYSIS 

In this section we utilize that the representation (|l.ip . (|1.2p of p„_i is of Cauchy- 
type which enables us to employ Riemann-Hilbert techniques. The key for this is 
that the function F„ defined in (jl.ip is the Cauchy transform of e""^ with respect to 
the contour 7. Hence, and this is the feature of the Cauchy transform that provides 
the link to Riemann-Hilbert problems, Fn is analytic in C \ 7 and the values of 
Fn{z) differ by e"'^'^'') as z approaches 567 from opposite sides. Moreover, these 
properties of together with the behavior of Fn{z) as z — s- cx) characterize F^ 
uniquely. More precisely, we will show in Lemma 13.11 below that Fn is the unique 
solution of the following scalar Riemann-Hilbert problem {RHP)i: 

Given an admissible contour^ and n G N. Seek an analytic function Y : C\7 ^ 
C such that 

(i) Y+{s) = r_(s) -t- e""^(") fors£-f , 

(ii) Y{z) for |z| -> 00 . 

Condition (i) is shorthand notation for the requirement that Y has continuous 
extensions from the interior of 7 (respectively from the exterior of 7) onto 7 which 
are denoted by 1+ (respectively Y-) and which satisfy relation (i). As mentioned 
above the question of existence and uniqueness of solutions for this Riemann-Hilbert 
problem is answered by the following 

Lemma 3.1. Fn as defined in is the unique solution of (RHP)i. 

Proof. Using (|1.2p it is easy to verify that Fn indeed solves {RHP)i. In order to 
prove uniqueness one first shows that the difference A of two solutions of (RHP) 1 
is continuous across 7 and hence entire. Liouville's theorem together with condition 
(ii) then implies A(2:) — for all z G C. □ 

Next we apply the method of constructing a local parametrix for the Riemann- 
Hilbert problem in order to resolve the difficulties at the critical point zq. We recall 
from Section [2] that the change of variables z t with z — X(it/y/n), z G J7o, maps 
the contour of steepest descent into the real axis and transforms n(\) to normal form, 
n(j){z) — —t^. This motivates the definitions 

(3.1) h{0 -L f ll^du, CeC\M, 

(3.2) Pn{z) h{-iV^X-\z)) , zgC/o\7- 

We have h+{t) = h^{t) + e^* for i G M, where h±(t) denotes lim,,^o+ h{t ± irj). 
This relation is a consequence of standard properties of the Cauchy transform. 
Alternatively, since h is the Cauchy transform of an analytic function e~" , one 
may also derive this relation using only the calculus of residues. Substituting s = 
\{it/y/n) we obtain 

{Pn) + {s) = {Pn)-{s) + e"*(^) for all s G f/o n 7 . 

Pn is thus a local solution of the Riemann-Hilbert problem {RHP)i in [/q- More- 
over, Pn is of a rather explicit nature, since h is related to the well studied comple- 
mentary error function [7] (see also (|4.5p below) and the Taylor coefficients of A can 
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be computed explicitly at ^ = to all orders from the defining relation (12. 3|) . Note, 
however, that this local solution Pn cannot be continued to the global solution be- 
cause X~^{z) has singularities outside of Uq. The procedure to take advantage of 
this local parametrix is strikingly simple. Choose e > such that the closed disc 
C/2e(l) is contained in Uq. We then set 

^ ^ ^ "'^"^-l r(z)-P„(z), forzeC/2e(l)\7 

where Y = F„ denotes the unique solution of {RHP)i. Observe that = 
TO_(s) for s g 7 n C/2e(l) since the jumps of Y and P„ cancel each other. We may 
therefore extend m within L'2e(l) to an analytic function. We denote this function 
with a slightly extended domain of definition by to. It is obvious that to again 
solves a Riemann-Hilbert problem. More precisely, denote 

Ti :-aC/2e(l), r2 :=7\C^2e(l), T Ti U . 

The dashed line in Figure [3] provides a sketch of F. It is straightforward to verify 
that the function to is a solution of {RHP)2: 

Seek an analytic function M : C \ F ^ C such that 

(i) M+is) = M_(s) - P„(s) for s e Fi \ F2 
M+(s) = M_(s) + e"''^('') for s G F2 

(ii) M{z) for |z| ^ 00 . 

As above M± denote the continuous extensions of M from the interior (+) and 
exterior (— ). We show in Remark 13.21 below that again we can use the Cauchy 
transform to solve the Riemann-Hilbert problem (RHP)2'. 

(3.4) TO(.)^^/-^(-)---(-)d.^^/ ^ds+± I ^^ds. 



2TTi 7p s — z 27ri Jp^ s — z 27ri 

The orientations of Fi and F2 are chosen to be counterclockwise so that the -|- 
side always lies to the left of the contour. Observe that all z G Ue{^) have at 
least distance e from F so that the contour of integration has been moved away 
from the singularity (s — z)~^ of the integrand. Moreover, the integral over F2 is 
exponentially small so that the asymptotic expansion of to is solely determined by 
the integral over Fi. Before we formulate our theorem on the asymptotics of Fn{z) 
for z near = 1 we present an elementary proof of (j3.4[) which does not make use of 
the fact that m is the unique solution of {RHP)2- Nevertheless, {RHP)2 provides 
a clear explanation why (|3.4p holds: The solution of this scalar Riemann-Hilbert 
problem is simply given by the Cauchy transform of the jump to-|_ — m_on F. 



Remark 3.2. One may show the first equality of p.4p by applying the calculus of 
residues to 

m(s) , f m(s) , f m(s) , 



where (Ti , a2 and as denote the closed curves shown in Figure [3] and by taking the 
limit as these curves approach F (dashed line). Note that condition (ii) in {RHP)2 
implies that m{s)/{s — z) has a vanishing residue at s = 00. 



We are now ready to state our main result. 
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Theorem 3.3. There exists e > and functions gj analytic in 11^(1) such that for 
any admissible curve 7 and any r € N we have 



(3.5) 



F„(z) = P„(z) + 




where the error term is uniform for z G [/^(l) \ 7 and the Taylor coefficients of all 
Qj are explicitly computable at zq — 1, e.g. 



9q{z) 



12 



(z-l) + 0((z-l)^) 



Proof. Choose e > such that the closed disc C/2e(l) is contained in UanU where 
Uq is defined above (|2.3p and U is determined by the curve 7 (cf. Definition 12 . iT ii) ) . 
From the discussion of the present section (see in particular p.3|) . (|3.4p ) it follows 
that the unique solution Fn — Y oi {RHP)i can be written for z G C^e(l) in the 
form 



F„(z) = P„(z) + m{z) = P,,{z) - / 

ZTTl Jp^ S — Z 



+ 0(e-"'=) . 



Observe that we may not use the calculus of residues to evaluate Jp ^s-'z since 
Pn is not a meromorphic function. Nevertheless, using 



,2r+l 



J=0 



C2'^+1(U-C) 



together with the fact that the Cauchy transform 2^ ^ — — ds is bounded on 
C\K (which can e.g. be proved using a contour deformation argument) one derives 
that for each r G N there exists a constant C such that 



Q2j 



< 



c 



ICI 



2r + l 
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for all C G C \ M. From the definition of P„ in terms of h we obtain 
(3-6) Pn{s) = ^ ( V + rest 





By choosing e smaller, if necessary, we can ensure that all Gj are analytic and 
non-zero in some open neighborhood of C/2£(l)- The term "rest" in p.6p is of order 
n^^ uniformly for s G dU2e{l). The calculus of residues then leads to the desired 
representation of Fn and the Taylor coefficients of gj can be derived from the Taylor 
coefficients of by explicit computation. □ 

Remark 3.4. Replacing C/2c(l) by ?7e/2(l) in (|3.3p and adapting the arguments 
above one obtains an expansion of F„(z) also for |z — 1| > e. Indeed, one can prove 
that there exist polynomials hj of degree 2j such that for all r G N we have 

(3.7) F„(z) = 

where the error term is uniform for e < |z — 1| < 2. Again the polynomials hj may 
be computed explicitly, e.g. hi{z) = —j^i^'^ + + 1). This result sharpens the 
statement of Theorem 12.21 It is also not difficult to see that the expansion p.7p 
holds outside shrinking circles n^" < \z ~ 1\ < 2 for any fixed < a < 1/2. In this 

case the error term has to be replaced by O (^ n'-\i-z\^'' ) ensure uniformity in z. 

4. Locating the zeros of the exponential sum 

The arguments and results put forward in this section will be presented in detail 
in a later publication, where we will construct for a specific admissible contour 7 
precisely n — 1 different solutions (together with their asymptotic expansions) of 
the equation 

(4.1) e"^(^) = Fn{z) 

in the interior of 7, provided n is sufficiently large. Since pn-i{nz) has at most 
n — 1 different zeros (jl.2|) implies that we have located all the zeros of pn-i- We 
will find the solutions of (|4.ip by showing for each 1 < k < n — 1 that the equation 

(4.2) G„(z):=^(z)-ilnF„(z)--^ 

n n 

has one solution in the unit disc Ui{Q). Here (p is defined as (p in (jl.ip with the 
only difference that the branch of the logarithm is now chosen in such a way that 
its imaginary part takes values in (0, 27r) rather than in (— 7r,7r). Existence and 
asymptotic expansions of the roots of (|4.2p are obtained by a standard procedure. 
First one constructs solutions n 

of 

, , , 27r?fc 

(4.3) An{z) = 

n 

for some approximation An of G„. Then a contraction mapping argument will be 
used to conclude that the original equation (|4.2p has a solution Zk,n close to ak,n- 
Since we can prove uniform bounds on the derivatives of the error terms in (|3.5p . 
()3.7p we obtain in addition an asymptotic expansion for Zk,n in terms of ak,n- 
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4.1. Zeros away from the critical point. In this case it suffices to use the crude 
approximation 

An{z) = 4>{z) = z - 1 - ln(z), 

where In denotes the branch of the logarithm described above. The solutions ak^i 
of (|4.3p all lie on the Szego curve Doo- One can show that the distance of ak^n and 
Zk,n is of the same order as An — Gn — -(ln_F!„) = 0((lnn)/n) (see Theorem 12. 2|) . 
Note that this result only holds for Zk^n that lie in a compact subset of C \ {!}. 
An extended version of Theorem 12.21 as discussed at the end of Section [5] allows to 
include also those solutions Zk,n that lie outside shrinking discs — 1| > n^" with 
< a < 1/2. The distance between ak,n and Zfc^„ is then of the order |2 - 
We state our result on the asymptotic expansion of 

Theorem 4.1. There exist polynomials Qj{x,y) of degree j in the variable y and 
of degree < 2j — 2 in the variable x such that for 0<P<l,n^<k< n/2, and 
r e N we have 



, , , ^ (afc,n,ln [V2^(afc,„ - 1)]) i 



In n \ / n 



n 



where the constant in the error term only depends on the choice of f3 and r. The 
polynomials Qj can be computed explicitly. For example, we have 

Qi{x,y) = -^y; Q2{x,y) = --^y^ + ]^xy - ^{x^ + lQx + 1). 

For r — \ (i.e. without correction terms) this result was first proved in [51 
(A.47)]. 

Faster rates of convergence can be achieved by using a better approximation of 
Fn in the definition of An that is provided by Theorem 12. 21 namely 

An{z) = (t){z) + - ln(V27m(l - z)) . 
n 

The faster rate of approximation comes at the price that the corresponding ap- 
proximate solutions afc,n now lie on n - dependent curves D„ rather than on the 
Szego curve. We leave the corresponding expansion of the zeros in terms of such 
approximate solutions ak,n for a later publication. Note that both [2l (A. 48)] and 
[1] also work with such better approximations. 

4.2. Zeros near the critical point. In order to formulate the result we first need 
to recall the definition of the complementary error function 

z oo 

erfc (z) = 1 -= / e^*-^ dt = —= / e"^^ dt for z g C, 

2 

where the path of integration of the latter integral is subject to the restriction arg(t) 
a with |a| < ^ as t — > oo along the path. It is well known (see e.g. [7]) that all 
the zeros of this function lie in the second and third quadrant of the complex plane 
(i.e. in the regions ^ < arg(z) < tt and — tt < arg(2;) < — f ) and that in the second 
quadrant there are countably many zeros of the complementary error function. We 
denote these zeros by Wfe, A: S N, and we can order them by modulus \wk\ < |?«fe+i|. 
Our result on the solutions Zk,n of (|4.2p reads as follows. 
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Theorem 4.2. There exist polynomials qj of degree j such that for < (3 < 1, 
1 < k < , and r ^ N we have 

_ -, , gj-(V2wfc) 




where the constant in the error term only depends on the choice of r and p. More- 
over, the polynomials qj may be computed explicitly, e.g. 

— 1 — Ix 

qi{x) = x; q2{x) = — - — ; q3(x) = — — — . 

Such a result was proved for r — 2 and /3 < 1/3 in [2, (A. 34)] sharpening and 
extending previous resuhs of [9], [4]. Since ~ for /c — > oo it foUows that 

the j-th term in the expansion above is of order (■^)-'^^- In particular, we have that 

\zk,n — 1| is of order ^J^■ Consequently, the expansion of Zk,n in terms of the zeros 
of the complementary error function holds for all zeros in shrinking circles of size 
j^(/3-i)/2_ Thus, for any e > and n^^^/'^^+^ < \zk.n - 1| < subsections 14.11 
and 14.21 provide different expansions for the solutions of (j4.2p . A short calculation 
shows that Theorem 14.21 yields better approximations in shrinking circles of size 
0{n^^/^), otherwise the approximation with ak,n on the Szego curve (see Theorem 
14. ip is more advantageous. 

We finish by explaining how the complementary error function enters the picture 
in the proof of Theorem 14.21 Introducing the auxiliary function 

v{0 := e^' erfc(C), 

one verifies that v{() — 2e''^ — v{—() holds for C g R and hence by the identity 
principle on all of the complex plane. Since v and the complementary error function 
have the same set of roots we obtain 2e'^''v{—Wk)~^ — 1. Elementary estimates on 
the arguments show that the correct value of the logarithm is given by 

(4.4) wl - ln{v{-Wk)/2) + 2'Kik = , for aU fc e N . 

Next we state the relation between v and the function h defined in p.ip . For all C, 
with positive real part the following relation holds 

(4.5) v{0 - 2h{zQ . 

To see this one verifies that g{C,) :~ 2h{i()e^^^ has the same derivative and the 
same limiting behavior for C — > c» as the complementary error function. Recall 
that the function h was used to define the parametrix P„ (|3.2p . Following the 
procedure described at the beginning of the present section and keeping the result 
of Theorem 13.31 in mind we approximate G„ by 

An{z) ■.^^{z)--\nPn{z). 

n 

Setting ak.n ■= ^{wk/\/n), using that cj) and (jj agree on the upper half plane, and 
that the real part of —Wk is positive, we obtain from ()2.3p . (14. Sp . and ()4.4p 

1 1 2Trik 

An{ak,n) = (tyiXi-Wk/Vn)) \nh{~iwk) = -[wl - \n{v{-Wk)/2)] = 

n n n 

satisfying (|4.3p as desired. 
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